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Abstract 

The author considers the planar rotational motion of the mathematical pendulum 
with its pivot oscillating both vertically and horizontally, so the trajectory of the 
pivot is an ellipse close to a circle. The analysis is based on the exact rotational 
solutions in the case of circular pivot trajectory and zero gravity. The conditions for 
existence and stability of such solutions are derived. Assuming that the amplitudes 
of excitations are not small while the pivot trajectory has small ellipticity the ap- 
proximate solutions are found both for high and small linear damping. Comparison 
between approximate and numerical solutions is made for different values of the 
damping parameter. 



1 Introduction 

Elliptically excited pendulum (EEP) is a mathematical pendulum in the ver- 
tical plane whose pivot oscillates not only vertically but also horizontally with 
7r/2 phase shift, so that the pivot has elliptical trajectory, see Fig. 1. EEP is 
a natural generalization of pendulum with vertically vibrating pivot that is 
one of the most studied classical systems with parametric excitation, so it is 
often referred to simply as parametric pendulum, see for example [1-4] and 
references therein. 

Dynamics of EEP has been studied numerically and analytically in [5]. Ap- 
proximate oscillatory and rotational solutions for EEP are the common exam- 

Email address: a_belyakov@inbox.ru 
(Anton O. Belyakov). 



Preprint submitted to Elsevier 



4 January 2011 



pies in literature [6-9] on asymptotic methods. Sometimes EEP is presented 
in a slightly more general model of unbalanced rotor [6-8], where the phase 
shift between vertical and horizontal oscillations of the pivot can differ from 
7r/2. EEP is also a special case of generally excited pendulum in [10]. 

The usual assumption for approximate solution in the literature is the small- 
ness of dimensionless damping and pivot oscillation amplitudes in the EEP's 
equation of motion. The author could find only one paper [11], where os- 
cillations of EEP with high damping and yet small relative excitation were 
studied. 

In the present paper we study rotations of EEP with not small excitation 
amplitudes and with both small and not small linear damping. Our analy- 
sis uses the exact solutions for EEP with the absence of gravity and with 
equal excitation amplitudes, when elliptical trajectory of the pivot becomes 
circular. ^ 

The paper is organized as follows. In Section 2 the dimensionless equation 
of EEP motion is derived. In Section 3 the exact rotational solutions and 

their stability conditions are obtained in the case, with no gravity and the 
circular trajectory of the pivot. In Section 4 first and second order approximate 
solutions are obtained by multiple scale method [13] for the close to circle 
trajectory of the pivot and high damping, where we assume that gravity is 
small or the frequency of excitation is high. In Section 5 for the same excitation 
and small damping second order approximate solutions are obtained with the 
use of averaging method [3,14]. In Section 6 both solutions in Sections 4 and 5 
are compared with the numerical solutions for different values of the damping 
parameter. 



2 Main relations 

Equation of EEP's motion can be derived with the use of angular momentum 
alteration theorem (e. g. [5]) 

where / is the distance between the pivot and the concentrated mass m; c is 
the viscous damping coefficient; 9 is the angle of the pendulum deviation from 
the vertical position; t is time; g is gravitational acceleration. 

^ When there is no gravity the model of EEP coincides with that of hula- hoop, see 
[12] and references therein. 
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Fig. 1. Scheme of the elhptically excited mathematical pendulum of length I. The 
pivot of the pendulum moves along the elliptic trajectory (dashed line) with semiaxis 
X and Y. 

It is assumed that the pivot of the pendulum moves according to the periodic 
law 



x — Xsm{Qt), y — Ycos{Qt), 



(2) 



where X,Y, and Q are the amplitudes and frequency of the excitation. 
We introduce the following dimensionless parameters and new time 



9 
V 



Y -X 



e — 



21 



Y + X ^ 



21 



nt, 



(3) 



in which equation (1) with substituted (2) in it takes the following form 

e + /3e + ii sin(T + e)^e sin(T - 6) - sin(^) , (4) 



were wc use formula Y cos(0 1) sin(^) + X sin(f2 1) cos{9) 



Y+X 



Y-X 



sin(nt + 9) 



sin(r2t — 9). Here the upper dot denotes differentiation with respect to 



new time r. 



3 Exact rotational solution when £ = and a; = 



Conditions £ = cu = mean that we find the mode of rotation for the circular 
excitation X — Y with absence of gravity = 0. In this case, we call equation 
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(4) the unperturbed equation 



(5) 



which has exact solutions 

e^eo-r, (6) 

where constants 9q are defined by the following equality 

sin(^o) = -, (7) 
A* 

provided that \/3\ < ji. 

To investigate the stabihty of these solutions we present the angle 9 as 9 — 
9o — T+r], where rj — r]{T) is a small addition, and substitute it in equation (5). 
Then linearizing in (5) and using equality (7), we obtain the linear equation 

fi + f3fi + ncos{9o)ri = 0, (8) 



According to the Lyapunov stability theorem, solution (6) is asymptotically 
stable according to the linear approximation if all eigenvalues of linearized 
equation (8) have negative real parts. Which happens when the following 
inequalities are satisfied 

P>0, iicos{9o)>0, (9) 



obtained from the Routh-Hurwitz conditions. Prom conditions (9), assump- 
tion /X > in (3), and equality (7), it follows for ^ > that the solutions 

9^9o-r, 9o^ arcsin + 2Trk (10) 



are asymptotically stable, while the solutions 

9^9o-r, ^0 = 7r-arcsinf^j +27rA; (11) 



are unstable, where k is any integer number. For negative damping, /3 < 
0, both these solutions are unstable. Prom now on we will assume that the 
following conditions are satisfied 

< ^ < (12) 
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which ensure the existence of stable rotational solution (10) as it is seen from 
(7) and (9). Indeed, in order to guarantee asymptotic stability /3 should be 
not only positive, but also strictly less than fi because of the second condition 
in (9), which can be transformed to inequality /xcos(6'o) = — > with 
the use of the positive root for /xcos(^o) from (7). 



4 Approximate rotational solutions when £ and co ~ ^/e 



We assume that values of e and u;^ are small of the same order of smallness, 
i.e. £ ~ a;^ <^ 1, so we can introduce new parameter w — ou^/e. 

One can deduct from (3) and current assumptions that either gravity g is 
small or the frequency of excitation Q is high with such damping c and mass 
m so that damping coefficient ^ ~ 1. 

All small terms are in the right-hand side of equation (4). To solve equation 
(4) we will use multiple scale method [13]. In this method general solution of 
equation (4) is assumed to be of the following form 

e = -T + eo + e9i + 8^92 + . . . (13) 



A series of time scales (independent variables). To, Ti, . . ., is introduced, where 
To = r, Ti = £r, .... So that 6' is a function of these time scales, 6'(To, Ti, . . .). 
Using the chain rule, the time derivatives become 



= T»o + £T»i + £^T>2 + . . . (14) 

or 

^ = Dl + 2eDoD, + ^^(DoA + Dj) + . . . (15) 

where D]^ — d^/dT^. Next the general solution and the time derivatives are 
substituted into equation (4), where sines are expended into the Taylor series 
with respect to e. By grouping together the terms with the same powers of e 
and equating to zero, a set of differential equations is obtained. 



DlOo + f3Do{eo - To) + ^ sin(^o) = 0, (16) 
Die, + (3Doei +fxcosi9o)9i = -{2DoD, + PDi){9o - Tq) 
+ sin(2To -eo)+w sin(To - ^^o), (17) 
Dle2 + (30002 + fxcos{eo)92 = fj.sm{9o)9l/2 - (2Do^i + /3^i)^i 

- {DoD2 + Dj + /3D2){9o-To) 

- (cos(2To -9o)+w cos(To - ^o)) ^i, (18) 
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where we denote w = cu^/s. We have aheady found solution (6) for equation 
(16) in the previous section. Here we consider the same stable regular rotations 
1:1 whose zero approximation is given by (10). Hence, Oq is a constant and 
consequently we have {2DoDi + (3Di){9o~To) = and (^0^2 + ^^? +/9^2)(^o- 
Tq) = 0. Thus, equations (17) and (18) can be written in the following way 

Die, + PDoe, + ^pt^ - n = sin(2ro - ^0) + sin(ro - ^0) (19) 

Dle^ + l3D^e2 + - 13^2 = (561/2 - {2D^D, + 6, 
- (cos(2To -9o)+w cos(To - ^0)) ^1, (20) 

where we denote ^sin(^o) = (5 and /icos(^o) = ^JW^^~W with the use of 
relation (7) and the second condition in (9). 

Jf..l First order approximation 

In consequence of conditions (12) non-homogeneous linear differential equation 
(19) can be presented in the following form 

= Ai cos(To) + Bi sin(To) 

+ A2 cos(2To) + B2 sin(2To) , (21) 

where A, = -w^/^i, B, = w^l - ^V/x^, A2 = B2 = ^1 - /^V/x^, 

lower index denotes harmonic number. Equation (21) has a unique periodic 
solution 

^i(To) = ai cos(To) + h sm(To) 

+ a2 cos(2ro) + 62 sin(2ro) , (22) 



Thus, the solution for (4) in the first approximation can be written as follows 



wherea,^-^^^V^;^-+^^\6,^ 
and 62 = z!^iE±(lV^^)£i. 

3/32+^2+4(4_2^^2_^2) 
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—e 



2 ^ cos(2ro - ^o) + (4 - Va*' - P') sin(2ro - ^0) 
3/32 + + 8(2 - Va*' - 



2/3cos(To-^o) + (l- V;^?^)sm(ro-^o) ^ ^ 

""^ /.2 + 1-2Va^2_^. ' (23) 

where constant Oq is defined in (10). 



4.2 Second order approximation 



Since (22) does not contain any constant of integration we set {2DqDi + 
(5Di)9i = in equation (20) and substitute in it an expression cos(2To — ^0) + 
■u;cos(To — ^^o) = Bi cos(To) — Ai sin(To) + B2 cos(2To) — A2 sin(2To) witli co- 
efficients defined in (21). Thus, equation (21) takes the following form 

Dle^ + + ^JiJ? - n = ^ + j2iK cos(nTo) + 5; sin(nro)) ,(24) 

^ n=l 



where coefficients in the right-hand side are the following 





(62^ + + 02^ + oi^j P + {Aibi + ^2^2 - ^202 - Biai) 


(25) 




(aia2 + 6162) P + (Aib2 + A2&1 - -6102 - B2ai)/2 


(26) 




(aj-bl) ^/2-(AA + Siai)/2 


(27) 


4 = 


{aia2 - 6162) P - (^2&i + Aib2 + Bia2 + B2ai)/2 


(28) 




(al - hi) (3/2 - (A262 + 5202)72 


(29) 


B[ = 


(ai^s - ^102) /3 - (Aifls - ^2^1 + Bib2 - B2bi)/2 


(30) 


B',= 


/3ai6i + {Aia-i - Bibi)/2 


(31) 


B's = 


{aib2 + ^102) P + (Aia2 + ^2^1 - -B1&2 - B2bi)/2 


(32) 




^a2&2 + (^2a2 -52&2)/2. 


(33) 



Periodic solution for equation (24) has the following form which is obtained 
from (59) in Appendix 



^2 (To) 



A' 



2V/^2 _ ^2 

^ (^2 _ V^^7^)a; + n(3B'^ 

n=l 



l)/32 + ^^ + n\n^ - 2V/^2 _ ^2) 



cos(nTo) 



E 



+ (n^ - v^^^^)i?; 

(n2 - l)/52 + + ^2(^2 _ 



=- sin(nro) , 



(34) 
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H=1, CO =0.3, e=0.2, p =0.5, Q =-0.34993, Z =-6.9004 
-0.7^ — 




Fig. 2. Angular velocities calculated from the first order approximate solution 
(23), second order approximate solution (35), and results of numerical simulation, 
when damping coefficient /3 is not small. 



where constant term is derived from (60) taking Aq = A'q/2. Thus, second 
order approximate solution can be shortly written in the following form 



g^-r + 0o + £0^{T)+e%{T), (35) 

where constant Oq is defined in (10), function 6i in (22), and function 62 in 
(34). In Fig. 2 it is shown how first and second order approximate solutions 
approach the numerical solution. 

In this section we have used the fact that all solutions in each time scale 

must converge to corresponding unique periodic solutions because damping 
/3 is not small. If it is not the case the multiple scale analysis becomes more 
complicated. In the next section we will tackle the problem of small damping 
(3 ~ y/e with the use of classical averaging technique. 
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5 Approximate rotational solutions when £ f« 0, a; ~ y^, and ^ ~ 

One can see in (3) that assumptions a; ~ ^ ~ y/i are valid for the high 
frequency of excitation f2 ~ with other parameters being of order 1. 

Another option is small gravity g ^ e along with small ratio c/m y^. 

After change of variable 9 — —t + equation (4) takes the following form 

+ sin(2T - v^i?) + \few sin(T - x/i-f?), (36) 

with small right-hand side, where we denote /3 = and as in the previous 

section w = u'^/e. With zero right-hand side equation (36) i? + fii} — (3 = 
would describe harmonic oscillations about vahic with frequency ^/JI. 
After Taylor's expansion of sines in the right-hand side of (36) about 1^ = 
we obtain the following equation 

ij + {fi + e cos(2r) + ew cos(r)) ^ - {3 
= -\/isin(2T) -I- •\/^wsin(r) — ^/e|3■^ -\- + ^(e), (37) 

which describes oscillator with both basic and parametric excitations. To solve 
equation (37) we will use averaging method [3,14]. For that purpose we will 
write (37) in the standard form of first order differential equations with small 
right-hand sides. First, we use Poincare variables q and ip defined via the 
following solution of generating system ^ + iid — 13 = Q which is (37) with 

£ = 

^ + gcos(V'), = -y///?sin(V'). (38) 



In Poincare variables equation (37) becomes a system of first order differential 
equations 



''"""^-/(r,?,^), (39) 



V^ = v7^-^/(r,g,V^), (40) 

where small function /(r, g, = \/i'/i(T, q, ip) + £/2(t, q, ip) + o{e) is the right 
hand side of (36), where 
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fi{T,q,ijj)^sm{2T) +wsm{T) + pq^/JIsm{'^jJ), (41) 
/2(r, q,i>) = - (cos(2r) + w cos(t)) + ^ cos(V') j 

3 

+ ^(^ + ?cos(V^)) , (42) 

meaning that /(r, g, t/^) = 0{^/e). Our next assumption is that ^/Jl — 1 ~ ^/£ 
which means that excitation frequency is close to the first resonant frequency 
of basic excitation component sin(T) and to the first resonant frequency of 
parametric excitation component cos(2r) in equation (37). Thus, system (39), 
(40) is transformed by = ^ + r to the standard form 



^ sin(C + T)/(T,?,C + T), (43) 



^ ^ ^ _ 1 _ J_ cos(C + r) /(r, C + r) , (44) 

with small right-hand side, where new slow variable C is often referred to as 
phase mismatch. 



5.1 First order approximation 



In the first approximation so called averaged equations can be obtained by 
averaging the system (43) , (44) over period 27r 



Q = / sm(Z + r) /i(r, g, Z + r) dr + o(v^), (45) 



2-K 



27r^ 





27r 



I cos(Z + r) /i(t, g, Z + t) dr + o{V~e). (46) 

where g and Z are the averaged variables corresponding to q and C,. After 
taking the integrals we have the following system 



g = cos(Z) - ^g + o(Vi), (47) 

Z = ^ - 1 + sin(Z) + o( Vi), (48) 

stationary solutions (g = 0, Z = 0) of which are the following 
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Z = arctan|| ^^^~ ^-" j +27rA: + o(l), (50) 

where we have substituted back w — io^jz and /3 = Symbol arctan 

stands for the principal value of the function on the interval from to tt. Note 
that the phase Z is determined to within 27r rather than tt, since the functions 
sin(Z) and cos(Z) obtained from equations (47) and (48) determine Z up to 
an additive term l-nk. Sohition of system (43-44) in the first approximation is 
g = Q + o(l), C = ^ + o(l) so the solution of (4) is the following 

^ = -T + - + ^flQ cos(Z + r) + o(x/i), 
A* 

which does not contain higher harmonics observed numerically. That is why 
we need to proceed to the second order approximation. 



5.2 Second order approximation 



In the second approximation averaged equations can be obtained as follows 



Q = f-Vicos(Z) + sin(Z)^ 



4 VA* J J 2VA^ 

+ ("4^ + 4^ '^^^^^^) ^ + ^^^^ 
Z = ^-l + ('^sin(Z) + ^ f- - 1^ cos(Z)^ 



8 i^v^^^V^vS""''^* 

stationary solutions (Q = 0, Z = 0) can be found numerically or with ab- 
sence of gravity (a; = 0) analytically. Solution of system (43-44) in the second 
approximation is the following 



q^Q + ^{- sin(T - ^) + I sin(2T + ^) + ^ sin(3T + Z)) 



2v^ 

+v^^ sin (2t + 2Z) + o(Vi), 



(53) 
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C = Z + (cos(t - ^) + I cos(2r + Z) + ^ cos(3t + Z) 

cos (2r + 2Z) + o{^/i), (54) 

Substitution of these expressions into (38) yields the second order approximate 
solution of (36) in the following form 



^?=^ + gcos(Z + r) + V£^sin(Z) (55) 
Vi^sin(r) _ v/^sin(2r) 

which after changes of variable 6 — — r + ^/e'& and parameters w — uj'^/e, 
(3 — f3/y/e results in the approximate solution of the original equation (4) 



e = -T + ^ + V^Q cos(Z + r) + sm{Z) 
fj, 4 

^u^_e^^ (57) 
4^ 3^ 

Agreement of solution (57) with the numerical experiment is shown in Fig. 3. 
We see that the amphtude of angular velocity oscillations is much higher than 
that for not small P in Fig. 2. 



6 Domains of applicability 



For asymptotic solutions in the previous two sections the parameter con- 
straints are more strict than those for the existence of stable exact solution in 
(12). Thus, for our analysis in Section 4 to be valid we must exclude cases when 
/3 — o(l) or/and V/^^ _ ^2 ^ ^(^ly ^j^g ^g^gg -^hgn p — ^(i) and fj, — 0{1) is 
studied in Section 4. Note that the case when both /3 — o(l) and // = o(l) 
meaning \/ ix^ — ^'^ — o{l) has already been studied in the literature, for ex- 
ample in the more general model of unbalanced rotor in [6]. The case when 
/3 = 0(1) and fi = o(l) is not feasible for asymptotic rotational solution. 
Indeed, with such assumptions the generating system 9 + /39 — has only 
constant solutions. These different cases are presented in the Table 1. 

To show quantitatively the limits of applicability of the assumptions in Sec- 
tions 4 and 5 we plot the absolute and relative angular velocity errors depend- 
ing on parameter (3 while excitation was constant A* = 1, see Fig. 4. 
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H=1, 03=0.3, e=0.2, P=0.01, Q =2.0348, Z =2.6838 



o numerical 
— approximate 




time X 



2*pi 



Fig. 3. Angular velocity 9 of the second order approximate solution (57) compared 
with the results of numerical simulatios in the case of small damping P 





smaD 


not; small B 


small ji 


studied in the literature 


no rotations 


not small ji 


studied in section 5 


studied in section 4 



Table 1 

Model assumptions on smallness of dimensionless damping /3 and dimensionless 
semiaxes half-sum of the ellipse along which the pivot of the pendulum moves in 
the problem to find pendulum rotations. In all cases we assume that dimensionless 
half-difference e of semiaxes is small as well as uP'. 



7 Conclusion 



The exact rotational solutions in the case of equal excitation amplitudes and 
zero gravity are obtained. The conditions for existence and stability of such 
solutions are derived. Based on these exact solutions the approximate solu- 
tions are found both for high and small linear damping, assuming that the 
amplitudes of excitations arc not small. Comparison between approximate 
and numerical solutions shows a good agrement for the damping values of the 
assumed order. 
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H=1, f0=0.3, £=0.2 



H=1, w=0.3, £=0.2 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

damping ^ damping p 



Fig. 4. Maximal absolute (left) and relative (right) deviation of analytically obtained 
angular velocities in (23), (35), and (57) from the numerically obtained 9. Solution 
(57) specially obtained for better approximation at small damping (3. Error in the 
right graph is calculated relative to the oscillation amplitude of angular velocity 9. 
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Appendix 



Non-uniform second order linear differential equation 

1^ + + yi^T^^ ^ An cos(nr) + 5„ sin (nr), (58) 



has solutions which obey the principle of superposition and contain only har- 
monics of the right-hand side in (58) 



=- sm(nr) (59) 



(n2 - l)/32 + + n^{n^ - 2^/^^ - 
where n = 0,l,2, Thus, solution with n = is the following 

401 = (60) 
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when n — 1 



(1 - VJ^iF^)A, + f^B 



+ 1 - 

-PA, + (1 - ^fiF^)B, 



- cos(r) 



sin(r) 



(61) 



when n — 2 



m 



(4 - ,/JF^)A2 + 
-2^^2 + (4 - y/JF^)B2 



cos(2r) 



sin(2T) 



(62) 
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